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Abstract 

Isotropic S = 1/2 quasi-one-dimensional antiferromagnets are considered 
within the bosonization method. The l/zj_-corrections to the interchain mean- 
field theory (where z± is the number of nearest neighbors in transverse to 
chain directions) are obtained for the ground-state sublattice magnetization 
So and Neel temperature T/v- The corrections to T/v make up about 25% 
of mean-field value, while those to Sq are small enough (especially in the 
three-dimensional case). The fluctuation corrections obtained improve con- 
siderably the agreement with the experimental data for magnetic-chain com- 
pounds KCUF3, Sr2Cu03 and Ca2Cu03. 
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I. INTRODUCTION 



Systems containing chains of magnetic atoms are investigated for a long time from both 
theoretical and experimental point of view. There exist many real compounds which are "al- 
most" one-dimensional (ID), i.e. have small interchain coupling. Here belong, e.g., KCuF 3 , 
Sr 2 Cu0 3 (spin S = 1/2), CsNiCl 3 (S = 1), CsVCl 3 (S = 3/2) etc. There are a num- 
ber of approaches which give a possibility to perform calculations for purely ID magnets 
(Bethe ansatz, exact numerical diagonalization, different versions of numerical renormal- 
ization group, quantum Monte-Carlo method etc.). At the same time, consideration of 
multichain problem with the use of these methods meets difficulties, so that theoretical ap- 
proaches are of interest which can adequately describe the situation in quasi-lD magnets in 
the presence of interlayer coupling and/or anisotropy. 

As for purely ID antiferromagnets, there is well-known theoretical result by Haldane |L| 
who mapped the spin-chain problem to nonlinear-sigma model (NLaM) and showed that 
the cases of integer and half- integer spins differ qualitatively (for a review see, e.g., Ref. 
0). For half- integer spins, the so-called topological ^-term in the effective action occurs 
which leads to unusual magnetic behavior of such chains. As follows from the Bethe ansatz 
solution for S — 1/2 (the same situation holds for any half- integer spin value), ground-state 
in this case already possesses quasi-long-range order. The excitation spectrum turns out 
to be gapless and spin correlators have a power-law behavior, but staggered magnetization 
is zero (the situation is reminiscent of the XY model below the Kosterlitz-Thouless point 
Tkt)- It is natural to suppose that in such a state the true long-range order is induced by 
an arbitrarily small interchain coupling J' and/or magnetic anisotropy. For the isotropic 
Heisenberg model, this problem was investigated within different theoretical methods. The 
interchain mean-field theory |3|-§] predicts for the ground-state staggered magnetization Sq 
and Neel temperature T N the results 

S otJ\jyj,T N <x\J'\ (1) 

and therefore indeed yields occurrence of long-range order at arbitrarily small \J'\. The 
behavior (0) contradicts to the standard spin-wave theory which does not distinguish between 
integer and half-integer spins and predicts in both the cases a finite critical value, J' c ~ Je~ nS , 
so that at | J'\ < J' c the quantity So vanishes and 

S (xln\J'/J' c \,T N (xSo^flFl (2) 

for | J' | greater but not too close to J' c . 

This contradiction was resolved within the renormalization-group (RG) approach 0J§-§] 
which showed that for inverse-length scales \x 3> J'J J the standard two-dimensional NLcrM 
scaling equations are applicable, and the spin-field scale factor indeed satisfies Z' 1 ! 2 oc 
In \i. At the same time, for half-integer spins at \x <C J 'J J one has Z~ l l 2 oc fx 1 / 2 . This 
means |||7]] that for both integer and half-integer spins and \J'\ ^> J' c we have the spin- 
wave behavior (||), while for half-integer spins and \J'\ <C J' c Eq.(|l|) holds. (We suppose 
here that for half-integer spins the renormalized coupling constant satisfies < g c where 
g c is the critical 3D coupling constant. Apparently, this inequality holds in the absence of 
dimerization, see Refs. H0|.) 
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In the extremely quantum case S = 1/2 we have J' c ~ J, so that | J'\ <C J' c in a broad 
region of \J'\. Therefore, one can conclude that interchain mean- field theory of Ref. @-[5| 
gives a qualitatively correct description of S = 1/2 quasi-lD magnets. At the same time, this 
theory does not take into account interchain fluctuations. In particular, the calculated value 
of the Neel temperature is not sensitive to space dimensionality of the system, although in the 
d = 1 + 1 case (both the dimensions are supposed to be spatial, but second one corresponds 
to the direction, transverse to the chain) we should have T/v = 0; for the d = 1 + 2 case the 
values of T/v turn out to be too high in comparison with experimental data. 

To obtain the corrections to interchain mean- field theory, we use the 1 / z_i_-expansion (z± 
is the number of nearest neighbors in directions transverse to the chain). This approach 
is similar to the expansion in 1/z (or inverse interaction radius 1/R), which has been used 
to improve the standard mean-field theory of Heisenberg magnets many years ago in Refs. 
PJlO[|. This approach is also equivalent to the spin-fluctuation approach in the theory of 
itinerant magnets by Moriya |[TT| . 

The plan of paper is as follows. In Sect. 2 we consider the bosonization of the system of 
interacting Heisenberg chains. In Sect. 3 we calculate fluctuation corrections to the interchain 
mean-field theory. In Sect. 3 we discuss the results and compare them with experimental data 
on magnetic chain compounds. In Appendix A the perturbation theory in J' is considered 
and the first-order l/zj_ correction to the mean-field value of Neel temperature are calculated. 
In Appendix B we demonstrate how the same results can be obtained more elegantly in spirit 
of the spin-fluctuation approach by Moriya. Finally, in Appendix C fluctuation corrections 
to the ground-state staggered magnetization are derived. 



II. THE MODEL AND ITS BOSONIZATION 

We consider the S = 1/2 isotropic Heisenberg model of quasi-lD antiferromagnet 

7~£ = J S^jSn+i^ + — J 1 ^ S ni jS n j (3) 

n,i n,<ij> 

where n numerates sites along the chains and i,j are indices of the chains, J > and J' are 
intra- and interchain exchange parameters respectively. We consider only the case \J'\ <C J. 
Each chain can be "bosonized" with the use of the standard relations (see, e.g., Ref. [O) 



S n , t = J l (x) + (-l) n n,(x) (4) 

where 

J*(x) = —d x (pi(x) 

jf(x) = —exp[±i/38i(x)] cos (3(pi(x) (5) 

7T 

are the cyclic vector current components and 

n z Ax) = — cos Pifi(x) 

TT 

nf(x) = — exp [±i(39i(x)} (6) 
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are their "staggered" analogs. Here A is the scale renormalization constant, <fi(x) is the 
boson operator, (3 = y/2%. 



Then we obtain the bosonized Hamiltonian in the form 



J'X 2 



J dx [ c os(A?i) cos((3(p i+s± ) + cos f3(9 i+s± - Bi)) (7) 



2tt 2 . , 

1,0 J_ 

where v = vrJ/2, Ilj is the momentum that is canonically cojugated to <fi, and Q{ satisfies 
d x 0i = —Ilj. The first line in (|?p corresponds to a system of separate chains and has the form 
of a standard sine-Gordon Hamiltonian. First term in (|7|) describes a free-boson system, and 
the second one arises because of Umklapp scattering of original fermions (which arises after 
applying the Jordan- Wigner transformation); this term is marginal and produces logarith- 
mic corrections to thermodynamic quantities P,P^TB|. Calculations (see Refs. PU*H}1) give 



g u /(2Tc) ~ 0.25. The second line of (^) describes the interaction between the chains. Note 
that only relevant terms are included in the second line since the marginal terms give smaller 
contribution (see Ref. 11131). 



III. MEAN-FIELD APPROXIMATION FOR BOSONIZED HAMILTONIAN AND 

1 /^-CORRECTIONS 

The simplest way of treating interchain exchange interactions is the mean-field approxi- 
mation Decoupling the interaction term 

cos((3(pi) cos(f3(p i+s± ) -> 2(cos(/3y? i+5± )) cos(^^) (8) 



we obtain 



where 



T~Cmf = t:Y1 f dx n ? + (dxPi) 2 ] + 9uYl J dxcos2/3ipi 

i i 

hM F y\ / dxcos(P(fi) (9) 

h MF = z ± J'\(cos((3(pi))/n, (10) 



zj_ is the number of nearest neighbors in the transverse (to chain) directions (z± = 4 for sim- 
ple cubic lattice). This approximation gives a possibility to reduce the multi-chain problem 
to a single-chain one in an effective staggered magnetic field. Introducing the function 

B(h ] T) = -(cos((3 Vl )) h (11) 

7T 

which should be calculated in the presence of the last term in ([]), we obtain the self- 
consistent equation for the sublattice magnetization S in the mean-field approximation in 
the form 
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Smf = B(z ± J'S M f;T) (12) 

Despite the Hamiltonian TImf © has an one-chain form, calculation of the function 
B(h;T) (which is an analog of the Brillouin function in the usual mean-field theory of 
Heisenberg magnets) at arbitrary T is a very complicated task. Scaling arguments suggest 
B(h; T) = h 1 / 3 f(h 2 / 3 /T) with some scaling function f(x). For g u = (in this case we have a 
standard sine-Gordon, or, equivalently, massive Thirring model) B(h;T) was calculated by 
Bethe ansatz in Ref . |I7| . However, in two following cases the calculation can be performed 
analytically: (i) T = where we have 0|| 

B(h; 0) ~ 0.677(V^) V3 [1 + ln(^/A)] 1/2 (13) 

where 

A ~ 2.085v 1/3 /i 2/3 

and (ii) h — > where 

B(h,T)=h X o{T). (14) 
Here Xo(T) is the staggered susceptibility of the system in the absence of h [§,0, 



where we have picked out the factor xo f° r the sake of convenience and 

AJ 1. , AJ „, J 1/2 



L(AJ/T) = C 



In— + - In In— + 0(1) 



(16) 



is the spin-field renormalization factor that arises because of the presence of the marginal 
operator; the single-chain numerical calculations [T3[] yield C ~ 0.15, A ~ 5.8. Thus one can 
see that the above-mentioned scaling function f(x) satisfies f(x) ~ x at x — > and /(oo) = 
const. The result ( ITj ) gives a possibility to calculate the value of T/v in the mean- field theory 
since for T — > T/v we just have /imf - > 0. Thus we obtain the equation |Q 

T^ = z ± J'x L(AJ/T*? F ) (17) 

We have included in ( |16|) a double-logarithmic term which was not taken into account in Ref. 
Q and modifies somewhat numerical results (see below). As discussed in the Introduction, 



the mean-field approximation (|17D is not quite satisfactory to describe experimental data. 
In particular, the values of Neel temperatures are considerably overestimated. 

The reason of this is that the mean-field approximation does not take into account the 
collective excitations which substantially contribute to the thermodynamic properties. Such 
excitations can be considered within the random-phase approximation (RPA). The RPA 
spin susceptibilities are given by [f|,[5| 

x fe ' w)= r^ro (18a) 

X ZZ (Qz, = J"^, V (18b) 



5 



where, for the square lattice in the direction transverse to chains, 



J'(Qx,q y ) = 2J'(cosq x + cosq y ), (19) 

Xo(q,uj) being the dynamical staggered susceptibility for the model (Jj|) and we have taken 
into account only staggered components of the susceptibility. Again, Xo(q,^>) is given by 
simple analytical expressions only in two cases: T = (for the results see Ref. M and also 



Appendix C), and h — > where we have [ 1Q , |T6 1 for both the susceptibilities 



Xo(q z ,uj) = f-) x (q z /T,u/T) 

v (k ,_i r(i/4 + tfc + )r(i/4 + ifc_) _u±k 

Xo(A; ' v) ~ 4r(3/4 + ifc + )r(3/4 + ifc_)' fc± - (20) 

As it should be, Xo(0, 0) = Xo(^)- 

Now we can calculate the spin-fluctuation corrections to interchain mean-field theory 
owing to collective modes. Similar to the case of the simplest mean-field approximation 
in the theory of Heisenberg magnets [PJnj, these corrections can be obtained within the 
1/z-expansion. Since we treat only transverse neighbors within the mean-field approach, 
one has to speak about the l/z^-expansion. To construct this expansion, we consider the 
perturbation theory in J'/ max(/i M ^,T) ~ l/z± (see Appendix A), which is an analog of 
expansion in J/ ma,x(hMF,T) ~ l/z for three-dimensional Heisenberg magnets [TO]. From 



this viewpoint, the above-discussed mean-field approximation is just the zeroth-order in 
l/zj_, so that the fluctuation corrections to this approximation can be obtained in a regular 
way. The leading (first order in l/z±) corrections come from the diagrams which include 
one RPA-interaction line. 

The details of calculations are discussed in Appendix A. For the Neel temperature we 
obtain to first order in l/z± 



T N = kJ'z ± xoL(A/T N ) (21) 



where 



°° 1 rl 1 

k= l\ '— I dr [ drV(r, t 



TV 2 



2xo 



-oo 



-F(r,r) + -G(r,r) 



oo 

rr f \ fdq z ^^ cosg^ + cosgy . 

Vlr, t) = / — > > — —z— exp(iq z r — 2mnr) (22) 

J 2n n gt^y 2 Xo ~ [cosq x + cosq y )xo{q z , 2mn) 



and F(r,r), G(t,t) are the four-point averages determined in Appendix A. The result (21]) 



differs from the mean-field result ([17]) by a factor of k which depends only on the lattice 
structure in the directions perpendicular to chains. Numerical calculation for d — 1 + 2 
case (simple cubic lattice) yields k ~ 0.70. Thus, with account of the function L(A/Tjv), the 
lowering of T N due to interchain fluctuation effects is about 25% of its mean-field value. For 
d — 1 + 2 the integral in ( p2"D is divergent and we have Tjv = 0. 

The same result (^TJ) can also be derived in a more elegant way within the spin-fluctuation 
approach by Moriya []lTJ (see Appendix B). 
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Corrections to the ground-state staggered magnetization are calculated in Appendix C. 
We have 



S = (0.677 - Q.311I)h]^ F (23) 

where the last term in the brackets represents the l/z± correction with 

_ / 0.038 d = 1 + 2 
1 ~ \ 0.193 d=l + l ■ [ ^ } 

Note that we do no take into account the logarithmic corrections owing to presence of the 
marginal operator, since there exists no simple ways of calculating dynamical staggered 
susceptibility at T = in the presence of such an operator. However, one can see that we 
have nearly 10% lowering of So for d — 1 + 1 and only 2% lowering for d — 1 + 2. Thus the 
fluctuation corrections to ground-state magnetization are much less important than those 
to the Neel temperature, and in the three-dimensional case they can be neglected. 



IV. COMPARISON WITH THE EXPERIMENTAL DATA AND CONCLUSION 

The results obtained enable one to perform quantitative comparison with experimental 
data on magnetic chain systems. Consider first the compound KCuF 3 with S — 1/2. Ac- 
cording to Ref. p(| , we have J = 406 K, S /S = 0.25. As discussed by Schulz [|]], this 
value of So corresponds to J' /J = 0.047, so that J' = 19.1 K. The simplest mean field 
approximation ( |T7| ) yields T/v = 47K. From fl2"T| ) we obtain T/v = 37.7 K which is somewhat 



lower in comparison with the experimental result of Ref. [20|, T/v = 39 K. Thus our ap- 
proximation slightly overestimates the effects of fluctuations, but improves reasonably the 
mean-field approximation. Contribution of the double-logarithmic term in (|T6|) makes up 
about 5 percents and improves the agreement with the experimental data. 

Another S = 1/2 chain compound that is widely discussed in recent publications is 
Sr 2 Cu03 which has the following parameters J = 2600K, T/v = 5K. Direct ex- 

perimental data for J' are absent, but using fl2~l"| ) and the experimental value of T/v we 
obtain J' = 1.85K. Then we have from (^) Sq/S = 0.042 which is in agreement with the 
experimental data (Sq/S < 0.05). 



For Ca2Cu03 the experimental parameters have following values pl|j22| : S = 1/2, J = 
2600K and T N = UK. From ([H]) we find J' = 4.3 K and S /S = 0.062. Taking into account 
above results for Sr 2 Cu0 3 we find that the latter value is again in excellent agreement with 



the experimental data p2f which give So(Ca 2 Cu03)/So(Sr 2 Cu03) = 1.5 ± 0.1. Thus the 



result (|21|) is sufficient to describe quantitatively real quasi- ID magnetic systems. 

In the isotropic quasi- ID magnets under consideration, the fluctuation corrections mod- 
ify only numerical factor in the expression for T/v- One can expect, however, that in the 
anisotropic case the form of functional dependence T^(J') will be also modified. The influ- 
ence of anisotropy on the Neel temperature will be considered elsewhere. Another interesting 
question concerns quasi-lD magnets with half-integer spins S > 1/2. As discussed in the 
Introduction, in this case there is a crossover from "usual" spin-wave behavior of staggered 
magnetization to non-spin-wave one. The expressions for T/v should be also changed because 
of this crossover. 
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Finally, despite the standard spin-wave theory yields a qualitatively correct description 
of integer-spin magnetic chains, the corresponding values of Neel temperatures are also over- 
estimated in comparison with experimental data. Thus calculation of fluctuation corrections 
for these magnets is also of interest. 
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VI. APPENDIX A. PERTURBATION THEORY IN J' AND THE DIAGRAM 
TECHNIQUE FOR SPIN OPERATORS 

In this Appendix we consider perturbation theory in J 1 for the field theory with the 
Lagrangian 

£=7^Y^J d 2 x(d(fi) 2 + g u Y^ J d 2 x cos 2/3^ - -h ^ J d 2 x cos((3(fi) 



J'\ 2 



J d 2 x [cos(/%) cos(f3(p i+s± ) + cos (3(9 i+s± - 0»)] (25) 



2tt 2 . . 

i,a± 

which corresponds to the Hamiltonian ([F|), the external staggered magnetic field h being 
introduced. In ([^D we have used the complex coordinate x = x + ivr. Farther in this 
Appendix we use the system of units where v — 1. Consider the calculation of staggered 
magnetization 

S = \(cos((3(pi)) /n 

The perturbation theory in J' is constructed in a standard way (see, e.g., Ref. p3| ). To 
obtain the series in J' we write down the expression in the path integral formalism 

-a A /£ycos(/3^(0))exp(-£[y]) 
7T fD<pexp(-£[- n 



To zeroth order in J' (i.e. at J' = 0) we have C = Co and 

S = B(h;T) (27) 
where the function B was introduced in (0). Expanding (p6| ) in J' we obtain 
^ = \f_D(pcos(J3(pi(p)) exp(-£ [y])(l - dm + C 2 nt /2 + ...) 



7T 



fD<pexp(-Co[<p])(l - L int + C 2 nt /2 



= -(COS(/%(0))(1 - C int + C 2 J2 + ...)) ,conn (28) 
7T 

where we have denoted (...) = / Dip.... exp(—C [{p})/ J Dip exp(— £qM) an d 
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n-1 



(cos(/3^(0))£^) ,c 



<cos(/?^(0))£L)o - E 



ran 



' m\(n — m)\ 



(cos(^(0))£-) <^ m )o 



(29) 



Each term in (|28|) can be represented by its own diagram; the diagram technique is the same 
as that for spin operators |p|,|T0|| (some elements of diagramm technique are shown on Fig.l). 
All diagrams are classified by powers of the parameter J'/ max(/i MF ,T) ~ l/z±. Diagrams 
of Fig. 2 have zeroth order in l/z±. The summation of these diagrams leads to a shift of 
the external magnetic field by the mean field: 



h — ► h = h + h 



MFi h>MF 



z ± J'S 



(30) 



(The same result could be obtained by eliminating the mean- field term directly in |25|). The 
diagrams of first order in l/z± (see Fig. 3a) have one RPA-interaction line (Fig. 3b). These 
are directly connected to the RPA susceptibilities ([18| ) by 



Thus we obtain 



V 



5 (q,w) = J'{q x ,q y ) + [J\q x , q y )fx + ' 2Z fe,w) 



J'{qx,q y ) 



l + 8 - J'(q x , q y )x ' zz {qz,uj) 



(31) 



where 



A 2 r 

Xo Z (Qz,v) = — J d 2 x (cos/3y?i(0) cos (3cpi(-x)}o,irexjp{-iq z x + iuj n T) 
Xt{q z ^) = -J rf 2 x(e^^(°)-^W]) exp(-^x + ^ n r), 
the two-operator irreducible average being given by 

(AB)i r — (AB) — (A)(B) 



(32) 



(33) 



and, following to [JTTJ] , we have introduced the correction 5 in the denominator to satisfy the 
self-consistency requirement. In the case T < T/v under consideration, this is determined by 
the condition [x + ~(0,0)] _1 = 0, i.e. 5 = z± J'xo ~(0, 0) — 1. Transforming (|3l|) back to the 
real space, 



V + -< zz (x) =TE E V + -> zz (q,zu n )exp(iq z x-tu n T) 



(34) 



Qx,q y 



we obtain for the sublattice magnetization (see the diagramms of Fig. 3a) 



S = B{h-T) + 



2f_ 

2tt 3 



d 2 xd 2 y [V zz (x - y)(cos/3</?i(0) cosf3ipi(x) cos /3<^(y)) ,i 



V + -(x - y)(cos/3^(0)e w(x) e- w(y) ) , i 



(35) 
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where 



(ABC) ir = (ABC) - (A)(BC) ir - (B)(AC) ir - (AB) ir (C) 



(36) 



(all averages are calculated with h — *■ h). Up to this moment, we did not use a concrete 
form of Co- As already pointed in the main text, the only case where the averages in (|35|) 
can be calculated analytically is the limit h — ► 0. In this limit we have 



A 



(cos/%(0) cos (3(pi(x) cos(3ipi(y)) 0iir 
h / (i 2 z(cos/3</?i(0) cos^(z) cos/3^(x) cos/5^(y)) ,ir 



(37) 



where 



(ABCD)i r = (ABCD) - (AD) ir (BC) ir - (BD) ir (AC) ir - (AB) ir (CD) l 



(38) 



and similar expression for transverse components; the averages in the right-hand side of 
are calculated at h — 0. Thus we have at h — 0, Kmf 



_ ^2 

S = ~r;h MF 
7T 



g? 2 z(cos P<fii(0) cos P<pi(z)) 



X 4 



2vr 4 



h MF / d xd yd z [V zz (x - y)(cos/%(0) cos/fy^x) cos/3</?i(y) cos (z)) ,-. 



+V + -(x - y)(cos/3^(0) cos^^(z) e ^«-^^) 0jir 
Note that the SU(2) invariance guarantees at h = 

d 2 x d 2 y d 2 z [(cos/5(^j(0) cos/3<^j(x) cos/3(^j(y) cos(3ipi(z)) 0t . 



(39) 



-3(cos/^(0) cos/3^(x)e i ^W- i ^W) 0lir 







(40) 



Calculating at /3 2 = 2n the averages in the right-hand side of (|3^) in the presence of the 
marginal operator g u cos2(3<pi (which produces logarithmic corrections) we obtain 



^ 1 



$ — —HmfL ( — 



A 



,2 1 1 , r2 f A 

d Z W)\ + ^ HmfL It 



x / d 2 xd 2 yd 2 zV^(x-y) 

k(y)?(z-x)| 



k(z)?(x-y)| 



+ 



k(zMz-y)| 



k( x MyM z - x M z - y) I k( z MyM z - x M x - y) I 

2 2 2 



k(z)?(x)?(z-y)?(x-y)| k(z)?(x-y)| |?(x)?(y-z)| |<?(y)<r(x - z)| 

+ h MF L 2 d) f d 2 xd 2 yd 2 zV + -(x-y)-—± 
4 T J kzkx 



|?(z)?(x-y)| 



xRe 



^(x)^(z - yfc(z - xMy) _ 1 
\| ?(x)?(z - y)c(z - x)q(y) 



(41) 



where 
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<j(x) = sinh(7rTx)/(7rT) 



(42) 



and L(A/T) is determined by fllEp. Introducing r = x — y instead of x and passing to the 
variables r = rT etc. we obtain the result 



2 



~F(r) + ±G(r) 



where 



7T 



Xo 



d 2 z T 



K z 



2.1184 



(43) 



(44) 



and 
F(i) 



d 2 y d 2 z 



+ 



\q(z)q(z-y)\ 



+ y)?(y)?(z - y - r)?(z - y)| |?(z)?(y)?(z - y - r)?(r)| 
?(y)fl(z-y-r)| 2 2 5 



G(r) 
In (E 



|?(z)s(r + y)?(z-y)c(r)| |?(z)?(r)| |<j(r + y)?(y - z)| k(y)?(r + y - z) | 
1 



<i 2 y <i 2 z- 
6f) we have used 



r Re 



g(r + yfc(z - y)q(z - y - r)c(y) _ 
^ + y)?(z - y)?(z - y - r)<f(y) 



<j(x) = sinh(7rx) 



(45) 
(46) 



(47) 



Finally, using the connection (RO) between the mean field and the staggered magnetization 
and collecting all corrections to the denominator analogously to the usual three-dimensional 
Heisenberg magnets 24 , we obtain the result (|2i"D of the main text. 



VII. APPENDIX B. SPIN-FLUCTUATION APPROACH TO 
THERMODYNAMICS OF QUASI-1D HEISENBERG MAGNETS. 

The results of previous Appendix can be obtained in a much more simple way with the 
use of the spin-fluctuation approach proposed by T. Moriya |nj for description of thermo- 



dynamics of itinerant magnets where the Stoner theory (which is an analog of the mean-field 
theory in Heisenberg magnets) turns out to be quite not satisfactory. To apply the spin- 
fluctuation approach, we represent Hamiltonian (]3|) as 

n = n + n int (48) 

^0 = J X! S n ,iSn+l,i — hMF y^A — ^) n+t S[ n 
n,i n,i 

1 <ij> 

where (AB) ir = AB — (A)(B). With the use of the Hellmann-Feynman theorem, we obtain 
for the free energy 
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(49) 



where Fo^mf) is the free energy corresponding to 7i - Using the RPA results (|T8|) one can 
find 

•T 7 = Fo(h MF ) +- J2 {in [l + J'(qx,q y )Xo~(qz,iu n ) + In [1 + J'(q x , q y )xo Z {qz, «*>„)]} (50) 



Differentiating with respect to /imf we readily obtain 



5 - S M f + - E 



J'(qx,g y ) 



dXo (qz,iu n ) 



1 + J'(qx,q y )Xo (qz,iuJ n ) 

J'(q x ,q y ) dxo z (q z ,iu r . 



dh 



1 + J'(q x ,q y )xo z (qz,iUn) 



dh 



(51) 



This is just the result (|35|) of Appendix A. Representing Xo(qz,iw n ) via boson variables, 
differentiating in h and calculating again the corresponding averages we return to (H3). 



VIII. APPENDIX C. GROUND-STATE FLUCTUATION CORRECTIONS IN 
THE ABSENCE OF MARGINAL OPERATOR. 

In this Appendix we consider ground-state corrections to the mean-field value of sublat- 
tice magnetization. We use the expression (|35|), or, equivalently, ( |51~D where, at T = (Refs. 



Xo 
Xo 



A\J'\uo 2 + v 2 q 2 + A 2 

Z'/Z A 2 

4|J'| w 2 + v 2 q 2 + 3 A 2 



(52) 
(53) 



with 



A ~ 6.175 \ J'\,Z'/Z 
S ~ 1.017| J'| 



0.49, 



(54) 



and huF — z^J'Sq (we neglect here the contribution of marginal operator). Differentiating 



(pl|) in huF (with account of implicit dependence of J 1 on h MF ) we obtain after some 
algebraic manipulations 



_ _ A OA 
S = S - ——I 



E 



4tt dh 

1 - r'/2) In — '— + (3 - Z'Y'/2Z) In 



/(3Z) 



(55) 



where we have used V = cosg for d = 1 + 1 and T r = (cosq x + cosq y )/2 for d = 1 + 2. 
Calculating the integral / numerically we obtain the result fl23|) of the main text. 
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Figure captions. 

Fig.l. Some elements of diagram technique for spin operators (for a detailed description 
see Ref. The first three irreducible averages are determined by ([33|), (pE|) and (|38|). 

Fig. 2. Diagramms for staggered magnetization to zeroth order in l/zj_ (mean-field ap- 
proximation). 

Fig. 3. (a) Diagrams of first order in l/z± for staggered magnetization (b) Equations for 
RPA interaction lines. 
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© = <S Z > (o~o) = <S i z S i z > ir 
( o •->•) = <Si Z Si + Si"> ir 
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